Structural transitions of ion strings in quantum potentials 
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We analyse the stability and dynamics of an ion chain confined inside a high-finesse optical res- 
onator. When the dipolar transition of the ions strongly couples to one cavity mode, the mechanical 
effects of light modify the chain properties close to a structural transition. We focus on the linear 
chain close to the zigzag instability and show that linear and zigzag arrays are bistable for certain 
strengths of the laser pumping the cavity. For these regimes the chain is cooled into one of the con- 
figurations by cavity-enhanced photon scattering. The excitations of these structures mix photonic 
and vibrational fluctuations, which can be entangled at steady state. These features are signalled 
by Fano-like resonances in the spectrum of light at the cavity output. 
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Crystals of singly-charged ions in traps are remark- 
able realisations of the phenomenon first predicted by 
Wigner [1 . The level of control experimentally achieved 
on these systems is impressive even at the quantum level 
and makes them promising candidates for several appli- 
cations ranging from metrology to quantum information 
processing [2-4 . Their versatility also offers the possibil- 
ity to study paradigmatic models of strongly-correlated 
many-body systems 1^,9^. 

Structural transitions in ion crystals have recently at- 
tracted renewed interest p!QHT4] . They are due to the 
interplay between the repulsive Coulomb interaction and 
the confining potential of Paul or Penning traps and can 
be controlled by varying, for instance, the aspect ratio of 
the trap potential. A prominent example is the linear- 
zigzag transition, that is classically described by Landau 
model [OT, while its quantum analogue belongs to the 
universality class of a one-dimensional ferromagnet [13] . 
Studies of quenches across the instability in the classi- 
cal regime showed that formation of defects follows the 
predictions of the Kibble- Zurek mechanism [151 [16]. 

The combination of ion and dipolar traps [171 |T8] opens 
further perspectives, such as the possibility of realising 
the Frenkel-Kontorova model [191 120] and of coupling 
ultracold atomic systems with ions [lH [22]. In Refs. 
[T71 [T8] the dipolar potential is classical, the quantum 
fluctuations of the electromagnetic field being very small. 
A very different scenario can be reached in presence of 
a cavity. Experiments with trapped ions in front of a 
mirror showed a mirror-mediated dipole-dipole interac- 
tion [23] and demonstrated the mechanical effect of the 
vacuum state of the electromagnetic field on a single ion 
[24] . The recent achievement of strong coupling between 
the optical transitions of ions forming a Wigner crystal 
and one mode of a high-finesse resonator [25l [26] sets the 
stage for the observation of novel self-organised struc- 
tures. In this regime, mechanical forces due to multiple 
scattering of a cavity photon can be infinitely long-ranged 
and may modify the structural stability even at the single 
photon level. The understanding of such dynamics can 
allow one to identify new control tools as well as to ac- 
cess new strongly-correlated states. The competition of 



long-range potentials of different nature, however, gives 
rise to a theoretical problem of considerable complexity. 




FIG. 1. The dipolar transitions of ions forming a chain couple 
with one mode of a high-flnesse optical cavity. The depth 
of the optical potential inside the resonator depends on the 
ions' transverse positions and on the strength of the laser 
pumping the cavity. This property gives rise to hysteresis 
in the structural conflguration and to quantum correlations 
between photonic and mechanical fluctuations. The scheme 
can be implemented in setups like the ones realised in [27l [28] . 



In this Letter we theoretically characterise structural 
properties of crystalline structures inside a standing- wave 
resonator, analysing in particular how the crystal struc- 
ture close to the linear-zigzag instability is modified in 
this environment. Figure [l] displays the main features of 
the system. A string of N ions of mass m and charge q is 
confined within an optical resonator by a radiofrequency 
trap, here described by a harmonic potential with axial 
and transverse frequencies Ua and cjt, respectively. The 
dipolar transition of the ions interacts with a mode of 
the cavity field which is pumped by a laser with strength 
T]. When cavity and pump are sufficiently out of reso- 
nance from the atomic dipole transition, the dynamical 
variables of ions and cavity field are described by the op- 
erators fj and pj , denoting the position and momentum 
of the center of mass of the j-th ion in the array, and 
by the annihilation and creation operators a and a^ of a 
cavity photon at frequency ujc. Their coherent dynamics 
are governed by Hamiltonian H = i^cav + ^ions + ^int- 
Here, Hc^v = —hAca^a — ihr]{a — a^) is the Hamilto- 
nian for the cavity mode in absence of atoms and in the 
reference frame rotating at the pump frequency ujp with 
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Ac = cjp — cjc, the Hamiltonian for the ions is given by 



,(1) 



and includes the kinetic energy, the trap potential Vtrap, 
and the Coloumb repulsion Fcoui- Finally, Hi^t = 
ha^ alio {ri^ ... ^vn) describes the interaction between 
ions and cavity field, with 
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Here, go is the strength of the coupling between the cav- 
ity and the ions' transition at frequency cjq, k the cavity 
wave vector, Aq = oOp — ooq the detuning of the pump 
from the dipolar transition, and a denotes the width 
of the cavity mode, that is generally smaller than the 
chain length: The number of ions coupling to the cav- 
ity mode, and hence contributing to [/q, is TVeff, with 
A/eff < N. Frequency Uq weights the nonlinear coupling 
between motion and cavity mode: It is the shift of the 
cavity frequency due to the ions inside the resonator, and 
conversely it is the mechanical potential exerted on these 
ions by a single cavity photon [29]. This term gives rise 
to mechanical effects that, for strong coupling, can be 
significant at the single-photon level. Incoherent effects 
arise from spontaneous decay of the dipolar transition at 
rate 7, cavity decay at rate 2/^, and thermalization of 
the ions' motion with an external reservoir which may 
be due to patch potentials at the trap electrodes [H [30] . 
We choose the detuning |Ao| to be the largest parameter, 
corresponding to the inequality |Ao| ^ 7,/^, |Ac|,^oV^7 
with n the mean intracavity photon number. In this 
regime the cavity-ion interaction is mainly dispersive and 
spontaneous emission can be neglected [31]. 

The Hamiltonian H formally agrees with the one de- 
rived for neutral atoms [29l[32- However, while in [291 132] 
the atomic interaction is a contact potential, here the ions 
repel via the long-range Coulomb repulsion. Therefore, 
in the first case the strength of the pump determines 
the quantum phase of the atoms in a non-trivial way 
[33^. For ions, on the other hand, quantum degeneracy 
is irrelevant but the strength of the cavity potential can 
substantially modify the crystalline structure. An effec- 
tive dispersive potential for the particles can be derived 
when retardation effects can be discarded. In this limit 
the cavity field is determined by the instantaneous set of 
positions of all ions coupled with the cavity mode and 
reads a = r]/{hz — iAeff), with Aeff = Ac — /7o, while the 
corresponding effective potential takes the form 



Kff = (^^^/^) arctan (— Aeff/Ai:) . 



(3) 



This potential gives rise to an effective long-range force 
between these ions whose strength scales with the coop- 
erativity C = ^o^eff/(^|Ao|). The ions' structure is then 
determined by the positions at which the total potential 
^tot = V;;rap + "^^Coui + "l^ff exhibits minima. 



Two situations can be identified depending on the 
value of the cooperativity C. For C <C 1, the potential in 
Eq. ^ approaches a classical potential whose depth is 
independent of the ions' positions. In this limit, when the 
ion string is parallel to the cavity axis (which corresponds 
to setting all values i/j = 0), the system provides a reali- 
sation of the Frenkel-Kontorova model with trapped ions 
[19 . When the string is instead orthogonal to the cavity- 
mode wave vector, as in Fig. [l] the optical potential gen- 
erates a transverse force. This force is symmetric about 
the chain axis when the chain is at a node or antinode 
of the cavity standing wave. Then, close to the linear- 
zigzag mechanical instability the optical potential shifts 
the critical value of the transverse trap frequency with re- 
spect to the free-space value ujtc [34 . In the following we 
shall assume that the equilibrium positions of the ions in 
the linear array are located at an antinode of the cavity 
standing wave. This can be realised, for instance, with 
the setups of Refs. [271 EH]- For blue-detuned pumps, 
with Ao > 0, the antinode is a maximum of the cavity 
potential and a mechanical instability thus appears at 
frequencies larger than ujtc^ while a red-detuned pump 
field will have the opposite effect [35]. This behaviour 
is significantly modified at large cooperativities, C > 1, 
where the cavity-mediated interaction between the ions 
becomes relevant. We consider Ac = 0, Aq > 0: In this 
case the intracavity field is minimum when the ions form 
a linear array, while it increases when their equilibrium 
positions arrange in a zigzag. This property can give rise 
to bistability of the linear and zigzag structures which 
can be observed in the mean value of the intensity /out 
of the field at the cavity output. An example of this be- 
haviour is shown in Fig. |2|a) where /out is plotted as a 
function of the pump intensity 77^ for a chain of 60 ions 
in a harmonic trap, assuming that the central region of 
the chain couples to the cavity mode and TVeff ~ 5.7. 

Further insight is gained by analysing the effective po- 
tential Vs of the zigzag mode which is the soft mode of the 
linear-zigzag transition in free space [10]. We first con- 
sider the simplest limit when the ions can be assumed 
to be equidistant (which describes the chain central re- 
gion or a chain in an anharmonic axial trap [36]) and the 
mode amplitude reads Xs = ^A — ^Yxj/^/N. Denoting 
bycjs 
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the frequency of the soft mode above 



the critical point in free space, the potential Vs when the 
ions are uniformly illuminated by the cavity field and for 
ujt > uJtc reads 
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The second term on the right-hand side of Eq. Q de- 
scribes the effect of the optical field. For r] = the 
cavity mode is in the vacuum state and the linear ar- 
ray is stable. The soft mode becomes unstable when 
the optical power is increased above the threshold value 
^th = N{l^C'^)/{4C)uj'^,K./ujR, with ujr = hk'^/{2m) the 
recoil frequency. For C > 1 parameters can be found in 
which both linear and zigzag configurations are stable. In 




FIG. 2. (a) Intensity /out at the cavity output for varying 
pump strength P = (?7/?7o)^, with r]o = ku)1/{2u)r). /out is in 
units of /i = Iout{P — 1) and is obtained for a configuration 
minimizing Vtot, numerically found using linear and zigzag 
chains as initial guesses, (b) Potential Vs as a function of 
the displacement x of the central ion (in units of l/Zc) for 
P = 130, 160, 190. The dashed lines indicate nodes at kx = 
di7r/2. The plots correspond to a chain of 60 ^°Ca^ ions with 
interparticle distance 4.3 /im in the central region, coupled to 
a cavity mode with wavelength 866 nm and transverse width 
<j — 14 /im (A^eff — 5.7). The other parameters are uJaj^i^ — 
0.1 MHz, cjt/27r = 2.26 MHz (the critical value is cjtc/27r :^ 
2.216 MHz), Ac = 0, Ao/27r = 500 MHz, kJ^i^ = 0.5 MHz, 
po/27r = 9.4 MHz, 7/27r = 10 MHz, resulting in C = 2. 



a finite chain the amplitude of the soft mode and the cou- 
pling of the ions to the cavity are not uniform along the 
chain. However, the potential energy as a function of the 
soft-mode amplitude gives similar qualitative results, as 
shown in Fig. [2|b). Here, for certain values of r\ the po- 
tential can exhibit three minima, corresponding to stable 
linear and zigzag arrays. We remark that the observed 
bistability is a consequence of the nonlinear dependence 
of the optical forces on the positions of the atoms within 
the standing- wave field. In the thermodynamic limit, if 
the region of the chain interacting with the cavity mode 
is finite, the effect of this coupling is a localized defect in 
the chain. For a finite system, nevertheless, forces acting 
on few ions can generate arrays close to zigzag configura- 
tions due to the long-range Coulomb repulsion [37[ [38] . 

The three metastable configurations can be observed 
when Ac = 0, Aq ^ 7 as the result of a cooling pro- 
cess due the strong coupling with the cavity [39]. In this 
regime the excitations of the emergent crystalline struc- 
ture reach a stationary state mixing photonic and vibra- 
tional modes. We analyse their behaviour by considering 
the coupled dynamics of the quantum fluctuations of field 
and motion. Be ^a = a — a the quantum fluctuations of 
the field about the mean value a, and {^Xj^^y^] the dis- 
placement of the ion localized at the equilibrium position 
f: determined by the balance of harmonic, Coulomb 
and mean optical forces. For convenience we introduce 
the normal modes of the crystal, that characterise the 
dynamics of the ions when the coupling with the quan- 
tum fluctuations of the cavity field can be neglected. Let 
K3 = En M^nBnoibn^bi) with C = ^, ^ and b^ (bl) the 



bosonic operator annihilating (creating) a phonon of the 
normal mode at frequency cj^, B^o = \/^/(27Tia;n), and 

M-^ the element of the orthogonal matrix relating the 
local coordinates with the normal modes. The dynamics 
of normal modes and field fluctuations are governed by 
the Heisenberg-Langevin equations [401 HI] • 

Set = (iAeff - f^)Sa - i ^ Cna(bn + b\^) + V^ain , (5) 

n 

bn = -{iuJn + ^n)bn - iCnip* 5a + a5a^) + \/2r^6in,n , 

(6) 

that include quantum noise on the cavity at rate n with 
corresponding input noise ain, and on the motion at rate 
Vn with input noise 6in,n, simulating the presence of a 
reservoir with which the ions' vibrations couple, such 
that (6L,n(^)^in,m(^0) = Nn5n,^S{t-t') [42^, with TV, the 
mean phonon number at the temperature of the reservoir. 
Vibrations and field fluctuations couple with strength 
c^a, where 
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and Qj = go cos{kxj)e~'^3 ^^'^^ \ The coefficients c, van- 
ish when all equilibrium positions are at field nodes, 
where gj =0. If the particles are located at antinodes, 
the coupling is determined by the derivatives in y direc- 
tion which are assumed to be much smaller than those 
along X {k(j ^ 1). Thus, for the chosen setup the cou- 
pling between vibrations and field fluctuations is stronger 
in the zigzag configuration, while it is a very small per- 
turbation when the ions form a linear chain. We remind 
that for the parameters considered the cavity field cools 
the normal modes coupled to it, so that cavity and crystal 
reach a stationary state [39 . 

We study the effect of this coupling in the spectrum 
at the cavity output, S{h') = (aout(^)^ <^out(^))/^0 7 with 
Iq = 2/^|ap the zero-order intensity of the output field 
and ttout = <^in + \f^a [40 . The steady-state spectrum 
exhibits negligible fluctuations in the linear phase, while 
in the zigzag configuration it reads 
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where the first term is the contribution due to the cou- 
pling of the quantum vacuum with the crystal vibrations, 
with Q(y) = '^n^n^n/[<^n + iln " ^^j^]? and the second 
is due to thermal noise coupling to the modes. The com- 
mon prefactor takes the form 



Soi'y) = 



K^ + {U + Aeff)2 



4Aeff^(iv)|a|' 



{k - ivy + A2g 



(9) 
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FIG. 3. Spectrum S{iy) of the field at the cavity output (in units of 5*0 = l/cJa) for a zigzag chain of 3 ions when only one 
ion at the chain edge, the right one in (c), couples significantly to the cavity mode (with Lamb-Dicke parameter ^ 0.1) and 
the ions' motion is cooled by the cavity field. The parameters are the same as in Fig. [2] expect for cua = /^ = 27r x 1 MHz, 
cjt = 27r X 1.57 MHz, and (a) C = 0.5, P = 1; (b) C = 3, P = 0.22; the mode width a is 0.65 times the interparticle distance in 
the linear array, so A^eff — 1.1. The equilibrium positions are the same in (a) and (b). The motional modes contributing to the 
spectral peaks are sketched in panel (c) (not to scale; the resonance corresponding to mode 2 is not visible, because this mode 
is too weakly coupled to the cavity). The Rayleigh peak is not shown. 



and gives a Lorentz curve when C <C 1. Its functional be- 
haviour is strongly modified when the cooperativity is in- 
creased: Then, motional and quantum noise do not sim- 
ply add up, but nonlinearly mix to determine the spectral 
properties of the output field. 

Figure [3] displays the spectra for a chain of three ions 
for different parameter choices: as C is increased the spec- 
tral lines change the relative heights, width, and shape. 
We first note the asymmetry in the spectra with respect 
to i/ = 0: This is due to the (weak) coupling of the ions' 
motion to the thermal bath [43 . The broadening at large 
cooperativity is a consequence of the vacuum input noise 
on the cavity field and indicates the rate at which the cav- 
ity cools the corresponding vibrational mode [44 . It is 
accompanied by the appearance of Fano-like resonances 
which result from the dispersive effect of the cavity back- 
action and are a signature of quantum interference in the 
fluctuations of motion and field [45 . This interference 
is due to quantum correlations established by the dy- 
namics described in Eqs. ([5])-(p|, which can generate 
entanglement between vibrational and photonic modes 
[46] . In fact, for the parameters of Fig. [s] (b) we find in 
the steady-state a logarithmic negativity of 0.15 [47 be- 
tween cavity and phononic excitations. We remark that 
the field at the cavity output allows one to monitor the 
stationary state in a non-destructive way, it can be mea- 
sured in existing experimental setups [48 and could be 
used to realise feedback on the ion crystal, for instance by 



means of an appropriate generalization of the procedure 
inRef. [49 . 

In summary, the structural properties and quantum 
fluctuations of an ion Coulomb crystal inside a resonator 
are strongly affected by the nonlinearity of the cavity cou- 
pling. This effect is particularly visible close to structural 
instabilities. We have focused on the linear chain close 
to the zigzag instability and shown that for large coop- 
erativity also the zigzag array can be made stable by the 
photon-mediated interaction between the ions (when a 
small region of the chain is coupled to the cavity, this cou- 
pling induces a localized defect in the chain with a zigzag 
form). This behaviour can be detected by hysteretical be- 
haviour of the intensity of the field at the cavity output 
as a function of the pump strength. The excitations of 
these structures reach a stationary state where phonons 
and photons are strongly correlated and can exhibit en- 
tanglement. At even larger cooperativity the dynamics 
studied in this work could be induced by one photon in- 
side the resonator, thereby providing an unprecedented 
control of many-body systems at the single-photon level. 
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